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Abstract 



A serious computational problem in the evaluation of hadronic collision pro- 
cesses is connected with the large number of partonic subprocesses included in the 
calculation. These are from the quark and gluon content of the initial hadrons, and 
' from CKM quark mixing. For example, there are 180 subprocesses which contribute 

to the VF+2jets process, and 292 subprocesses in VF+3jets production at the LHC, 
even when quarks from only the first two generations are taken into account. 

We propose a simple modification of the rules for evaluation of cross sections 
and distributions, which avoids multiplication of channels from the mixture of quark 
states. The method is based on a unitary rotation of down quarks, thus, transport- 
ing the mixing matrix elements from vertices of Feynman diagrams to the parton 
distribution functions (PDF). As a result, one can calculate cross sections with sig- 
nificantly fewer subprocesses. For the example mentioned above, with the new rules, 
one need evaluate only 21 and 33 subprocesses respectively. The matrix elements 
of the subprocesses are calculated without quark mixing but with a modified PDF 
convolution which depends on the quark mixing angle, and on the topologies of 
gauge invariant classes of diagrams. The proposed method has been incorporated 
into the CompHEP program and checked with various examples. 



Introduction 

The evaluation of cross sections for hadron collisions typically involves a large number 
of subprocesses. The reason for this is the quark and gluon content of the initial hadrons 
and the Cabbibo-Kobayashi-Maskava (CKM) mixing of the down quarks. At Tevatron 
and LHC energies in many cases 5 flavors, u, d, c, s and b, give a sizable contribution 
through the corresponding parton densities. Note that subprocesses sometimes contribute 
only because of the non-diagonality of the CKM matrix. An example is the subprocess 
ud — > sdW~^ which contributes to the W + 2 jets production. Furthermore, the number of 
subprocesses increases with the number of quarks in the final state, in particular, because 
the final state jets produced by various subprocesses are indistinguishable Q The process 
of VF-boson and jet production (an important background to various Standard Model 
and "new physics" processes) exemplifies the problem: 180 subprocesses contribute to 
W + 2jets and 292 subprocesses toW + Sjets production if only the quarks of first two 
generations and only QCD diagrams are taken into account. 

The huge multiplication of channels stands as a real computation problem despite the 
fact that the matrix elements of some subprocesses may have similar analytical structures. 
Indeed, one should separately organize Monte Carlo integration and/or event generation 
for each subprocess separately because the convolution with the parton distribution func- 
tions (PDF) is flavour dependent. This is a present-day problem for automatic calculations 
of collision processes (see e.g. and references therein for review of this new computation 
approach) . 

In this letter we propose a new method which simplifies the flavour combinatorics and 
reduces the multiplication of channels due to the mixture of quark states. Implementation 
of the proposed algorithm in Monte Carlo codes is straightforward. It has been incorpo- 
rated into the CompHEP program ||^ and checked for many examples. In the method 
the quarks of the first two generations are taken to be massless and do not mix with the 
third generation which is obviously a good approximation for many processes in the Teva- 
tron and LHC energy range. The method is based on a rotation of down quarks, thus, 
transporting the mixing matrix elements from vertices of subprocess Feynman diagrams 
to the parton distribution functions. The complete set of subprocess diagrams are divided 
into gauge invariant classes, and depending on the topology of the class, new computation 
rules are formulated. The method drastically simplifies the calculation of matrix elements 
but leads to a modification of the PDF convolution procedure. We demonstrate the power 
of the proposed technique with the example of W + 2jets production at LHC. Here only 
21 subprocesses have to be evaluated with the new technique, which should be compared 
with 180 subprocesses in the standard approach. 

The paper is organized as follows. In Section 1 we discuss CKM diagonalization at 
the level of Feynman amplitudes and introduce basic notations. Here the topologies of 
the gauge invariant classes of diagrams are defined, for which in Section 2 new PDF 
convolution rules are formulated. The W + 2jets production process at LHC is discussed 
in Section 3. In the Appendix we briefly comment on the implementation of the new rules 
in the CompHEP program. 

^An important exclusion is a 6-quark jet which can be distinguished from Hght quark jets with a 
micro- vertex detector. 
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1. Feynman amplitudes 

Let us consider a parton subprocess with a quark in the initial state. There are two 
possible topologies of any Feynman diagram with the corresponding quark line: 1) the 
quark line goes through the diagram from the initial to final state - "scattering topology" , 
and 2) the quark line connects both initial state partons - "annihilation topology". In 
Fig.[^ and Fig.^ one can see examples of scattering topology, and diagrams of annihilation 
topology are shown in Fig.Hb,c. Here the interior of the circle represents the diagram, 
and the solid line represents a quark line (quark current). An important point is that 
according to the theorem proved in the two sets of diagrams with the scattering and 
annihilation topologies form gauge invariant classes with respect to the SM gauge group. 
Thus, one can change the computational rules for each class independently. Note, that 
a specific parton subprocess could contain only scattering topology diagrams, or only 
annihilation topology diagrams, or both. For the process ud —>■ ddW~^ shown in Fig.^ 
both annihilation and scattering diagrams contribute. For our further consideration it is 
necessary to determine whether the quark current is charged (CC) or neutral (NC). This 
is done by counting the quark vertices with VF-bosons. The charged current involves an 
odd and the neutral current an even number of W vertices. Both cases might include an 
arbitrary number of vertices with gluons, photons and Z-bosons. 

Two approximations are used in our method. We neglect: l)the masses of quarks from 
the two first generations, and 2) the mixing with the third generation. Thus, the 3x3 
CKM mixing matrix reduces to the 2x2 Cabbibo matrix: 



where "i^c is the Cabbibo angle. 

The zero quark mass limit is a good approximation in many applications, in particular, 
when one evaluates matrix elements of hard subprocesses. Indeed, the energy scale, e.g. 
the partonic collision energy, for Tevatron-LHC hard subprocesses is (9(100) GeV and 
higher. Then, minimal energy-like cuts on jets (momentum transfer, invariant masses 
of jet pairs etc.) are of order (9(10) GeV or more. This means that the momenta in 
the propagators are of the same order and mass effects can be neglected in matrix ele- 
ments. At the same time, the flavour dependence of parton distributions is sizable and 
can substantially affect observables. 

The second approximation (no mixing with the third generation) also works well in 
many applications. The Vtb mixing matrix element is very close to unity and the non- 
diagonal CKM matrix elements for bottom and top quarks are very small. Of course, these 
elements are responsible for important physical phenomena, such as B-meson oscillations 
and rare decays of heavy mesons. In these cases the proposed method can, obviously, not 
be applied. Another reason why the second approximation is reasonable is that t and b 
quarks produce flnal state objects in a detector which are very different from light quark 
jets: the b quark jet has a secondary vertex, and the t-quark appears as a heavy narrow 
resonance. 

The starting point of our consideration is the diagonalization of the quark mixing 
matrix in vertices. This means that rotated down quark states are used in Feynman 
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rules, d! = d cos 'dc + s sin 'dc and s' = —d sin 'dc + s cos 'dc rather than d and s states 
being eigenstates of the mass matrix. It is worth recalhng that all electroweak vertices 
are diagonal over the isodoublets ( ^/ ) and ( ) • For example, the charged current 

(CC) electroweak vertex is diagonal in this rotated basis, 

< ■ Jlic = K ■ u\ Y V.jd,'^ = ■ u\ 7^ 4^ . 

Here rf^ = Vijdj (where = 1,2 are the generation indices, so d = di and 5 = ^2)- As 
a result, elements of the CKM mixing matrix do not enter in the matrix element if one 
calculates in terms of these rotated down quarks. 

Let us consider an electroweak model with the only one generation of quarks (referred 
to as the EWud model), and denote generalized up and down quarks by Qu and respec- 
tively. Then, Feynman amplitudes in the EW model with two generations and Cabbibo 
mixing can be evaluated in the EWud model with a single quark generation. Indeed, in 
the EWud model, the CC vertex represents the Lagrangian term W^lqu^'j^qd'" , and its 
contribution to the amplitude is analytically the same as the contribution of the standard 
CC vertex (remember that we have neglected quark masses). The only difference is a 
multiplication of some amplitudes by mixing matrix elements. There are only two generic 
variants for the case of scattering topology diagrams illustrated in Fig.[T|: CC in Fig.|I|a 
and NC in Fig.|l]b. One can easily see that the rotation of down quarks results in only 
one factor Vij (or Vj~^ in conjugated cases) in front of the each CC line, while the mixing 
matrix elements cancel each other inside the NC line due to the unitarity of the mixing 
matrix J2k ^ik^ ^kj = ^ij- One should stress that the summation over internal flavours in 
propagators is performed when the amplitude is evaluated. In the Fig.|l| only one W ver- 
tex is shown in the CC case and two W vertices in the NC case. However the statements 
are correct for the general CC and NC cases. Obviously the same conclusions are valid in 
the case of CC and NC diagrams of the annihilation topology as shown in Fig.pD,c. 

An important difference between the annihilation and scattering topologies appears at 
the level of squared diagrams when the convolution with the parton distribution functions 
is performed. 

2. Squared diagrams and PDF convolution 

The next step is to square the matrix element and convolute with the PDF's. Here 
we derive four basic rules for the technique under discussion. 

Scattering topology. 1st Rule 

Let us consider Feynman diagrams with the scattering topology. The second m-state 
can be a quark, an anti-quark or a gluon. We denote the sum of this type of Feynman 
diagrams as T>sc- An example is presented in Fig. ^ for the case of a CC upper quark line 
and an NC lower quark line. After squaring of diagrams of this type and convoluting with 
PDF's one can write the contribution to the cross section in the form (for the example in 
the Fig. |a): 

l^^scP =^ j dXldX2fdX^l) fuj{x2) Wijkn? = 

ijkn 
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dxidx2 



■ ik 



where by M and M. we denote matrix elements evaluated in the standard EW theory 
and in the EWud model respectively. Due to the unitarity of the mixing matrix one has 
VikVf^^ = Sii, and the first rule can be written as 



dxidx2 [fd{Xl) + fsixi)] [fu{x2) + fc{x2)] \M\ 



(1) 



This means that one can evaluate squared diagrams of the scattering topology class 
with only one quark generation, but should convolute the corresponding gauge invari- 
ant squared matrix element with modified structure function(s) - a sum of down (or up) 
PDF's. 



Annihilation CC topology. 2nd Rule 

In the case of the annihilation topology the CC and NC cases lead to different rules for 
convolution with the PDF. We start from the CC annihilation case which is generically 
shown in Fig. ^d. Recall that this variant occurs only if the quark line has an odd number 
of W vertices. We denote the sum of this class of Feynman diagrams as 'D'^^ . When one 
squares this sum of diagrams and convolutes with the PDF, one gets: 



CC\2 



j dxidx2 fdX^i) fuj 



X2 



dxidx2 



j:Uixi)fn,{x2)V,,V^^' 



l-MI 



In contrast to the scattering case here one can not use the unitarity condition to cancel 
two elements of the mixing matrix because the summations over the indices, i and j, also 
include the structure functions fd^ and fuy Thus, the second rule can be written in the 
following form (see the example in Fig. [^d): 



w 



CC\2 



dXidX2 [ fd{xi) fu{x2) COS^ ^^c + fs{xi) fc{x2) COS^ + (2) 

+ fd{xi) fc{x2) sin^'d, + fsixi) fu{x2) sin^^^J \M\^ . 



where we have explicitly expressed the mixing matrix elements in terms of the Cabbibo 
mixing angle. 

One can see that the annihilation-type contribution to the squared matrix element is 
convoluted with non-factorizable products of PDF's. 



Annihilation NC topology. 3d Rule 

We denote the sum of all diagrams with the annihilation NC quark line as "Pf*^*". The 
generic example is shown in Fig. In this case the quark line has an even number of 
M^-boson vertices. In the same way as above one can derive the following formula: 

\'^a''? =^ Y J dXldX2 fdXxi) fd,{x2) = 
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= J dxidx2 

Therefore, the third rule can be written in the form (if the m-quarks are down): 

\1^a^\' =^ I dXrdX2 [Uxi) h{x2) + fs{Xl) f-s{x2)] \M\^ • (3) 

with the obvious generahzation for up quarks. 

Interference of Vsc and Va topologies 

In the general case diagrams of both, scattering and annihilation, topologies could con- 
tribute. The interference between gauge invariant diagram classes with these topologies is 
also gauge invariant and, therefore, can be independently convoluted with the PDF. As in 
the above derivations the PDF convolution for the interference of Vgc diagrams with 
diagrams is given by the same formula as for (2nd Rule), and for the interference 

of T>sc with T)'^^ by the corresponding formula for \T>^^\'^ (3rd Rule). 

The final state quark-antiquark line. 4th Rule 

Finally let us consider Feynman diagrams where a quark line connects two out-states, 
as in Fig. §d. Here the summation over generation indices does not involve parton distri- 
bution functions and can, therefore, be performed explicitly. The result is that the con- 
tribution of the corresponding squared diagrams, evaluated in the EWud model, should 
be multiplied by two for both the CC and NC cases. Indeed, according to Fig. |^d each of 
these summations give in the squared diagram the factors: Y^ij ^ij ^ji = 2 in the NC and 
J2ij Vij Vji^ = Y.ij ^ij = 2 in the CC cases. 

Note that the 4th Rule is valid not only in cases where the quark loop in squared 
diagrams connects out-state(s) and never passes through 2n-state(s), but also for each 
quark loop in next-to-leading corrections. 

Of course, the 4th Rule is valid only under the assumption that the fragmentation 
of the four light quarks and antiquarks leads to indistinguishable jets. If one includes 
nontrivial fragmentation functions, e.g. for a c-quark, all the above rules have to be 
modified. We do not present here the corresponding formulas which however could be 
easily derived. 

3. Test: W + 2jets production at LHC 

In this section we illustrate the proposed technique with the example of + 2jets 
production at LHC. Here a total of 180 subprocesses contribute in the standard technique 
when all four light quarks contribute separately and the CKM matrix is present in W- 
boson vertices. With the new technique only 21 subprocesses need to be evaluated. 

We shall not discuss this example in full detail but present the results of a numerical 
test for the ud ddW~^ subprocess with permutations of quarks within pairs {u,c) and 
{d,s). Note that in this example three rules are used: 1st, 2nd and 4th. In the standard 
technique 12 subprocesses are involved and the corresponding contributions to the cross 
section are collected in Tab. |^ for two values of the kinematical cut on the transverse 
momenta of the final partons and the W boson, pj^^ > p9p and > with pj, = 20 



\M\ 
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or 200 GeV. The cross sections were calculated in the Standard Model with the averaged 
values for CKM matrix elements and quark masses taken from the Particle Data Group 
P . For calculations the CompHEP program has been used, and the accumulated MC 
error in all cases was less than 0.6%. We neglect the contributions of subleading diagrams 
with electroweak boson propagators, calculating the cross section to leading aa^ order. 

Using the new method one should evaluates only one subprocess in the one quark 
doublet EWud model, Quqd QdQdW^ . The corresponding Feynman diagrams are shown 
in Fig. ^. This subprocess is of the mixed type where two gauge invariant classes of 
diagrams: the annihilation CC topology (Fig. ^) and the scattering topology (Fig. ^), 
contribute. The 1st Rule is used for the squared scattering topology contribution. The 2nd 
Rule is used for the squared annihilation CC topology contribution with a multiplication 
by the factor 2 according to the 4th Rule. Finally, the 2nd Rule is used to evaluate 
the interference between the two gauge invariant classes of Feynman diagrams. We have 
calculated all three contributions to the cross section of the subprocess quQd QdQdW^ 
using CompHEP code in which the EWud model and Rules 1-3 have been implemented (see 
the Appendix for details). The results (see Tab. |^) for the total rate, a{pT > 20GeV) = 
112.95 pb and a{pT > 200GeV) = 0.30213 pb, are in an agreement with the "standard" 
calculations of Tab. within the statistical error of less than 0.6%. 

Conclusions 

We have shown that hard collision processes at hadron colliders can be evaluated in 
an economical way, greatly reducing the number of contributing subprocesses. The pro- 
posed computational technique can be applied only if the quark masses of the first two 
generations and mixing with the 3rd generation can be neglected. These assumptions 
are valid for most applications at the Tevatron and LHC. In the proposed technique the 
Standard Model with a single generation of up and down quarks is used and a squared 
matrix element is evaluated without involving elements of the mixing matrix. The re- 
sulting matrix squared matrix element is convoluted with modified parton distribution 
functions according to formulas given above as Rules 1-4. Each of these Rules corresponds 
to a gauge invariant class of squared diagrams. 

In Rule 1 the squared matrix element is convoluted with fu{x) + fdx) or fd{x) + fs{x) 
for an m-state of up-type or down-type quarks respectively However, in the cases of Rules 
2 and 3 the squared matrix element is convoluted over Bjorken variables Xi and X2 with 
a non-factorizable function. When applying the 2nd Rule the function 

[fd{Xl) fu{x2) + fs{Xl) /e(x2)] COS^ + [fdixi) /g(x2) + fs{Xi) /^(Xz)] slu^ , 

or a similar function where the substitution d d and u ^ u are made, is used. In the 
case of the 3rd Rule the convolution is with either 

fu{,Xi) fu{x2) + fc{,Xi) fc{x2) or, 
fdixi) f^{x2) + fsixi) fs{x2) , 

depending on which quark-antiquark pair occurs in the m-state. 
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Appendix 

The Rules derived in this letter have been implemented in the CompHEP code v. 33. 
New model, referred to above as EWud, was created on the base of the SM where only one 
quark generation, up and down quarks denoted as Qu and qd was kept without any CKM 
matrix elements. The masses of these generalized quarks were set to zero. Then, the 
option for numerical convolution of squared diagrams with parton distribution functions 
was modified in accordance with Rules 1-3. The code of this version of CompHEP is 
available from the Web 0. 

In the most general case the user has to subdivide the whole set of squared diagrams 
into two parts: 1) l^^scP, and 2) jX'aP plus the interference diagrams 2Re(Vsc-1^a)- Each 
of these parts should be calculated separately. In particular, for each part the user has to 
set the variable "PDFfactor" in the menu option "User menu" as follows: PDF factor = 1 
for l^scP (referred to in CompHEP as the "t-channel" case), and PDFfactor = for 
\Va\^ + 2Re{Vsc ■ VI) (referred to in CompHEP as the "s-channel" case). The program 
automatically recognizes which Rule, 2nd or 3rd, should be used in the latter case. 

Rule 4 (multiplication by a factor 2) has to be applied by hand if the corresponding 
quark line is presented in a squared diagram. 

One should note that this is not a completely automatic realization of Rules 1-4. This 
is planned for future CompHEP development. 
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Figure 1: Generic scattering topology diagrams with: a) a single W vertex and b) two W 
vertices on the quark line. No summation over the generation indices i and j is implied. 
Here and in the next figure the external lines on the left and right sides of the diagram 
correspond to in- and ou^-states respectively. 
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a) 




Figure 2: Generic diagrams with different topologies in wiiicli tlie quarlc lines connect in- 
and ou^-states. Here i and j are the generation indices of the m-quarks, while k and n 
are the generation indices for out-quarks. Diagram d) includes a quark line connecting 
two out-states. 
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^^Qd Qd < ' < Qd 



Figure 3: Feynman diagrams for the subprocess quQd QdQdW^ which contribute to 
+ 2jets production at LHC. Case a) corresponds to the V^'-^ topology of diagrams, 
while case b) corresponds to the Vgc topology. 



9 



Tables 



Subprocess 



cross section, pb 



= 20 GeV 



= 200 GeV 



ud 
ud 
ud ■ 
us 
us ■ 
us ■ 
cd ■ 
cd - 
cd ■ 
cs ■ 
cs - 
cs - 



■ ssW+ 
dsW+ 
ddW+ 

■ ssW+ 
dsW+ 
ddW+ 
ssW+ 
dsW+ 
ddW+ 
ssW+ 
dsW+ 
ddW+ 



9.368 

2.315 

55.84 

2.113 

32.67 

0.3748 

0.0739 

4.601 

0.3206 

4.093 

0.1528 

0.9842 



0.01746 

0.007807 

0.1733 

0.005301 

0.09074 

0.0005516 

0.00003557 

0.004845 

0.0002915 

0.003090 

0.0001359 

0.0003853 



In total 



112.90 



0.30391 



Table 1: Contributions, in pb, of different qq' channels to the + 2jets production 
cross section at LHC In addition to cuts on transverse momenta, p^^* > pj^ and pj^ > p^^, 
the following cuts were applied: on pseudorapidity Irjjetl < 5 and on the jet separation 
in the r] — ip plane R{jet, jet') > 0.5. The statistical accuracy of the MC calculations is 
better than 0.6%. 



Classes of squared diagrams, 
used Rules 


contribution to cross section, pb 


p^ = 20 GeV = 200 GeV 


iP.cP, 1st Rule 
Ica/D^^P, 2nd + 4th Rules 
2Re{V^^ -V*^), 2nd Rule 


89.76 0.2621 
21.62 0.03692 
1.549 0.003111 


In total 


112.93 0.30213 



Table 2: Contribution, in pb, of the subprocess quqd (ldQdW~^ to the + 2jets 
production cross section at LHC, evaluated with the help of Rules 1,2 and 4 using the 
of EWud model. The same cuts were applied as in Tab. The accuracy of the MC 
calculations is better than 0.6%. 
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